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The energy conditions and the Dolgov-Kawasaki criterion in generalized f{R) gravity with ar- 
bitrary coupling between matter and geometry are derived in this paper, which are quite general 
and can degenerate to the well-known energy conditions in GR and f{R) gravity with non-minimal 
coupling and non-coupling as special cases. In order to get some insight on the meaning of these 
' energy conditions and the Dolgov- Kawasaki criterion, we apply them to a class of models in the 
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^1 1. Introduction 

^jQ. Recent astrophysical obscrvations[l, 2] have indicated that the expansion of our Universe is accelerating at the 

present time. In principle, this phenomenon can be explained by either dark energy (see, for instance, Ref. [3] for 
reviews), in which the reason of this phenomenon is due to an exotic component with large negative pressure, or 
modified theories of gravity [4]. Unfortunately, up to now a satisfactory answer to the question that what dark energy 
is and where it came from has not yet to be obtained. Alternative to dark energy, modified theories of gravity is 
' extremely attractive because the cosmic speed-up can be easily explained by the fact that some sub-dominant terms, 
like may become essential at small curvature. Under some additional conditions, the early-time inflation and 

, late-time acceleration can be unified by different role of gravitational terms relevant at small and at large curvature. 
f{R) gravity is one of the competitive candidates in modified theories of gravity (see, for instance, Refs.[5, 6] for 
reviews). Here f{R) is an arbitrary function of the Ricci scalar R. One can add any form of R in it, such as l/i?[7], 
lni?[8], positive and negative powers of i?[9], Gauss-Bonnet invariant [10], etc. It worth stressing that considering some 
additional conditions, the early-time inflation and late-time acceleration can be unified by different role of gravitational 
terms relevant at small and at large curvature. However, /(i?) gravity is not perfect because of containing a number 
of instabilities. For instance, the theory with 1/R may develop the instability [22]. But by adding a term of R^ to 
this specific f{R) model, one can remove this instability[8, 9]. For more general forms of f{R), the stability condition 
/" > can be used to test /(i?) gravity models [25]. 

This paper is organized as follows. In section 2, we give some fundamental elements of generalized f{R) gravity 
models with arbitrary matter-geometry coupling. In section 3, the well-known energy conditions, namely, the strong 
energy condition (SEC), the null energy condition (NEC), the weak energy condition (WEC) and the dominant energy 
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condition (DEC), in the generalized f{R) gravity models, will be derived. In order to get some insight on the meaning 
of these energy conditions, we apply them to a class of models. Furthermore we rewritten them in terms of parameters 
of the deceleration (q), the jerk {j) and the snap (s) and then use the rewritten WEC to restrict a special f{R) model. 
The instability of generalized f{R) gravity models with arbitrary matter- geometry coupling will be studied in section 
4. Last section contains our summary. 



2. GENERALIZED f{R) GRAVITY MODELS WITH ARBITRARY MATTER-GEOMETRY COUPLING 

A more general model of f{R) gravity, in which the coupling style between matter and geometry is arbitrary and 
the Lagrangian density of matter only appears in coupling term, has been proposed in Ref. [11]. Its starting action is 

I + G{L,n)f2{R)]V^d^x, (1) 

where we have chosen k = 8ttG = c = 1. fi{R) {i = 1,2) and G{L„i) are arbitrary functions of the Ricci scalar R 
and the Lagrangian density of matter respectively. When f2{R) — 1 and G{Lm) ~ Lm, we obtain the general form of 
f{R) gravity with non-coupling between matter and geometry. Furthermore, by setting fi{R) — R, action (1) can be 
reduced to the standard General Relativity (GR). 

Varying the action (1) with respect to the metric yields the field equations 

Fi(i?)i?^, - yi{R)g^u + (5p.n - V^V,)^^i(i?) = -2G{L,n)F2{R)R^, 

-2{g^,D - V^\/,)G{L„-,)F2{R) - h{R)[K{L^)L„, - G{Lm,)]g^,, (2) 
+/2(i?)/v(i„)T^„ 

where □ = g^'^V^V^, Fi{R) = dfi{R)/dR (i = 1, 2) and K{L,n) = dG{L„i)/dL,n respectively. The energy-momentum 
tensor of matter is defined as: 

In this class of models, the energy-momentum tensor of matter is generally not conserved due to the appearance of 
an extra force [12]. 



3. ENERGY CONDITIONS IN THE GENERALIZED f{R) GRAVITY MODELS WITH ARBITRARY 

COUPLING BETWEEN MATTER AND GEOMETRY 



3.1. The Raychaudhuri Equation 



Many models of f{R) gravity have been proposed, which can be restricted by imposing the so-called energy 
conditions [13]. These energy conditions were used in different contexts to derive general results that hold for a 
variety of situations. Under these energy conditions, one allows not only to establish gravity which remains attrac- 
tive, but also to keep the demands that the energy density is positive and cannot fiow faster than light. Below, we 
simply review the Raychaudhuri equation which is the physical origin of the NEC and the SEC [14]. 
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In the case of a congruence of timelike geodesies defined by the vector field u^^, the Raychaudhuri equation is given 
by 

^ = -1:0^- ^^.^o""" + ^M-^^" - (4) 
ctr o 

where i?^^, 6, a^u and lo^^i, are the Rieci tensor, the expansion parameter, the shear and the rotation associated with 
the congruence respectively. While in the case of a congruence of null geodesies defined by the vector field fc*^, the 
Raychaudhuri equation is given by 

do ^ 

From above expressions, it is clear that the Raychaudhuri equation is purely geometric and independent of the 
gravity theory. In order to constrain the energy-momentum tensor by the Raychaudhuri equation, one can use the 
Rieci tensor from the field equations of gravity to make a connection. Namely, through the combination of the field 
equations of gravity and the Raychaudhuri equation, one can obtain physical conditions for the energy-momentum 
tensor. Since = a^^a^'^ > (the shear is a spatial tensor) and uj^ij ~ (hypersurface orthogonal congruence), 
from Eqs. (4) and (5), the conditions for gravity to remain attractive {dO/dr < 0) are 

R^^u'^u'' > SEC, (6) 

R^yk^'k" > NEC. (7) 

Thus by means of the relationship (6) and Einstein's equation, one obtains 

R^.uu^u'' = [T^, - ^gf,,)u^u'' > 0, (8) 

where T^i, is the energy-momentum tensor and T is its trace. If one considers a perfect fluid with energy density p 
and pressure p, 

Tfj,u = {p + p)U^Uy-pg^^, (9) 
the relationship (8) turns into the well-known SEC of Einstein's theory, i.e., 

p-f3p>0. (10) 
Similarly, by using the relationship (7) and Einstein's equation, one has 

T^.fc^r > 0. (11) 
Then considering Eq.(9), the familiar NEC of general relativity can be reproduced as: 

p + p>Q. (12) 

3.2. Energy conditions 

The Einstein tensor resulting from the field equations (2) is 

G^v — R^i/ — ■zQfj.uR ~ T^l'^ , (13) 
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where the effective energy-momentum tensor Tl^l^ is defined as follows: 

T^l^ - 7pi^{i.9^.[/i-(/{ + 2G/^)i?]-(.9^,n-V^V,)/( ^^^^ 

where fi — fi{R) [i = 1, 2), G ~ G(L,„) and the prime denotes differentiation with respect to the Ricci scalar R and 
the Lagrangian density i,„ respectively. Contracting the above equation, we have 

T'" = 7r+W{2[/i - ifi + 2G/^)i?] - 3n/( ^^^^ 

-eOG/^ - 4/2(G'i™ - G) + /aG'T}, 
where T — g^'^T^^^. Thus, wc can write in terms of an effective stress-energy tensor and its trace, i.e., 

R^.u=T;lf -]^g^,T^ff. (16) 

In order to keep gravity attractive, besides the expressions (6) and (7), the following additional condition should 
be required 

Note that this condition is independent of the ones derived from the Raychaudhuri equation (i.e., the expressions (6) 
and (7)), and only relates to an effective gravitational coupling. 
The FRW metric is chosen as: 

ds'^ = dt^ ~ a'^{t)dsl, (18) 

where a(t) is the scale factor and contains the spacial part of the metric. Using this metric, we can obtain 
R — -~6{2H^ + H), where H — a{t)/a{t) is the Hubble expansion parameter, and = a{t)d{t)5 (/i, v ^ Q) are 
the components of the affine connection. 

By using the relationship (6) and Eq. (16), the SEC can be given as: 

mfu^^u^ - ^T^ff > 0, (19) 

where we have used the condition g^j^^u^u'^ = 1. Taking the energy-momentum tensor T^i^ to be a perfect fluid (i.e., 
Eq.(9)) and considering the condition (17), we obtain 

P + 3p^j^[h- {f{ + 2G/^)i?] + 3^{HR + ii) 
+f^'R^G + f^RG) + 6j§jr{G'L^f^ + f^'RG) 



h 

-J,(G'i„-G)>0, 



where the dot denotes differentiation with respect to cosmic time. This is the SEC in f{R) gravity with arbitrary 
coupling between matter and geometry. 

The NEC in f{R) gravity with arbitrary coupling between matter and geometry can be expressed as: 

Tl^l^kt^k" > 0. (21) 
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By the same method as the SEC, the above relationship can be changed into 

L„,G'f^ + 2,f!{RG'K, + f^'R^G + f^RG)- (22) 
^(G'L;„/^ + /^'i?G) >0. 

From above discussions, it is worth stressing that by taking G{Lm) = Lm and rescahng the function f2{R) as 
1 + A/2(i?) in expressions (20) and (22), we can obtain the SEC and the NEC in f{R) gravity with non-minimal 
coupling between matter and geometry, which are just the results given in Ref. [15]. While by setting f2{R) — 1 and 
G{L„i) ~ Lm, we can derive the SEC and the NEC in f{R) gravity with non-coupling, which are just the same as 
the ones in Ref. [16]. Furthermore, when fi{R) — R, the SEC and the NEC in general relativity, i.e., p + 3p > and 
p + p > 0, can be reproduced. 

Note that the above expressions of the SEC and the NEC are directly derived from Raychaudhuri equation. However, 
equivalent results can be obtained by taking the transformations p — > p'''^^ and p — > p'^^-^ into p + 3p > and p+p > 0. 
Thus by extending this approach to p — p > and p > 0, we will give the DEC and the WEC in f{R) gravity with 
arbitrary coupling between matter and geometry in the following. 

By means of Eqs. (14) and (18), the effective energy density and the effective pressure can be derived as follows: 

P"" = 7F5G7I^5[/i - (/I + 2G/^)i?] - SHRfi' - m{G'Lmf^+ ^^^^ 
fl,'RG)-h{G'L,n'G) + f2G'p}, 



P"" = jr^^{~^[fi-if{+2Gn)R] + {2HR + R)fi' + ,fi"R^ 

+2{G"Ljf^ + LmG'f^ + 2/^'7?G'l;„ + f;{'R^G + /^'i?G) (24) 
iH{G'L„J^ + f^'RG) + h{G'Lm - G) + /sG'p}. 

Then, the corresponding DEC and WEC in f{R) gravity with arbitrary coupling can be respectively written as: 

P-P+j^[h- {.f'l + 2Gf^)R] - {5HR + R)^- 

J^^"^ ~ /2G" (^"-^rn /2 + -^mG" /2 + 2/2' -RG'i„i + 

f^'R^G + f^'RG) - ^{G'Lmf^ + f^RG) - ^ 
{G'Lra - G) > 0, 



{G'Lmf^ + f^RG) - ^{G'Lm - G) > 0. 

We show that by taking G{Lm) = Lm and rescaling the function f2{R) as 1 + A/2(i?), above expressions are the 
DEC and the WEC in /(i?) gravity with non-minimal coupling between matter and geometry, which arc just the 
same as the ones in Ref.[15]. While by setting f2{R) ~ 1 and G{Lm) — Lm, the results given by us are the DEC and 
the WEC in f{R) gravity with non-coupling, which are consistent with the results given in Ref. [16]. Furthermore, 
when /i(i?) = R, the DEC and the WEC in general relativity, i.e., p — p> and p > 0, can be reproduced. 
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3.3. Energy Conditions for a Class of Models 

In order to get some insight on the meaning of the above energy conditions, we consider a specific type of models 
where fi{R) and f2{R) arc taken as 

h{R) = R + eR-, ^2^^ 
/2(i?) = ai?™, 

In the FRW cosmology, the energy conditions can be written as 

^^^i^{-[G(L„Oa„ + A] I R +aj > (28) 

where A^B and Cm.„ depend on the energy condition under study and we take e = (— l)"e and a = (— l)"a due to 
the fact that for a FRW metric one has R <{). For the SEC, one finds 

^SEC ^ G'{Lm)[L,n + 3mR-\Lln + HL,n) + 6mLmRR-^m-l)] ^ ^ 

■ 2 (29a) 
+3mL„ G"{L„^)R-\ 



B^^^ ^-{p + Zp), (29b) 

Cf^^ = (?i- l)[3i?7ii?-2 + l + 3i77ii?-2i? + 37ii?-3i?2(„_2)]. (29c) 
For the NEC, one obtains 



A^^^ = [mL„,G'{L„,)R-^ - HmLmG'iLm)R-^ 

+2mL„,RG'{L^)R-^{m - 1)] + mL^G" {L,n)R'^^ , 



^DEC ^ G'{L,n)[-L,n-mCnR^^ -5HmL„^R-^ 

+ (1 - m)2TnL,nRR-^] - mL„i^G"{L,n)R-\ 



(30a) 



B^^^^-ip + p), (30b) 
C^sc ^ {n-l)n[RR-^ + HR-^R + R-^R^{n-2)]. (30c) 

For the DEC, one has 



(31a) 



BDEC (31b) 



C^^^c ^ (l-n)[i?ni?-2 + m_5^^^-2^_^„^-3^2(^_2)]. 



(31c) 
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Finally, for the WEC, one gets 

j^WEC ^ -L,n-6HmL^R-\ (32a) 
BWEC ^ (32b) 



^WEC ^ {l^n){^ + 3HnR-^R). (32c) 

Given these definitions, the study of all the energy conditions can be performed by satisfying the inequality (28). 
Note that all the energy conditions depend on the geometrical parameters. It means that for different models, the 
energy conditions can be satisfied by choosing them properly. 

For models given by Eq.(27), the condition for keeping gravity attractive(GA), i.e. inequality (17), also can be 
obtained from inequality (28) by taking 

^""^ = 2k I ^ I""' (33a) 



= 0, (33b) 

C,f^ = nR~\ (33c) 

This means that inequality (28) also stands for the condition that ensures gravity remains attractive for models given 
by Eq.(27). 

In the following, we use energy conditions to restrict a special f{R) model also in the FRW cosmology. The Ricci 
scalar R and its derivatives can be expressed by the parameters of the deceleration (q), the jerk (j) and the snap 
(s)[17], namely, 

R = -6H^{1- q), (34a) 
R^-6H^j -q-2), (34b) 
R= -6H^{s + q^ + 8q + 6), (34c) 



where 



1 a 1 a 1 a' 

'7 = -77J-, j = 77J-' s = — — . (35 



Thus, the energy conditions (20), (22), (25) and (26) can be rewritten as: 



P + Sp-jk^[fl+ QH^ifi + 2G/^)(1 - q)] - l^J^H\j + s + q^ + 7q + A) 
+108^i/6(j - g - 2)2 + q^[G"lJj^ + L^G' - Uf^'H^J - q 
-2)G'L,n + 36f!,"H^{j - q - 2fG - &f^H\s + q^ + 8q + 6)G] + 6^ 
[G'l;„/^ - Qf!lH\j - q - 2)G] + MG'L„, - G) > 0, {SEC) 



(36a) 
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p + p- 6HHj +s + q^ + 7q + 4)^ + 36^H%j -q-2f + j^ 
[G"Ljf^ + L^G'f^ - 12/^'i73(j - q - 2)G'Lm + 36f^"H^{j - q 
-2)2G - 6/^'i74(s + q^ + 8q + 6)G] - ^[G'Lrnf^ - Gf^H^j ~ q 
-2)G] > 0, (NEC) 

P-P+J^[.fi+ ^H\f[ + 2G/^)(1 - q)] + 6H\5j + s + q^ + 3q 
-4)^ - 36^H%j -q-2f- j^[G"Ljf^ + L^G' - Uf^H^ 
(j-q- 2)G'L„ + 36/f i/^C? -q- 2fG - 6f^'H^is + q^+8q + 6)G] 
-^[G'L„J^ - Gf^H^ij -q-2)G]- ^{G'L„, - G) > 0, {DEC) 

2^[/i + + 2G/^)(1 - 9)] + - 9 - 2)^ - &Hj^ 

[G'L,nf^ - Gf^H^j - q - 2)G] - ^(G'L„, - G) > 0. (WEC) 

To exemplify how to use these energy conditions to constrain the /(i?) theories of gravity, we consider a special 
model with ,fi{R) = R, f2{R) = oiR^ and G{L,n) = Lm = — p[18]. Since there has been no reliable measurement for 
the snap parameter (s) up to now, we only focus on the WEC. Under the requirement f'{R) > for all R and taking 
Hq = 70.5[19], the WEC (36d) in this particular case is 

0.3Bn2-0.3n(l + B) + l>0, (37) 

where B ^ {j - q - 2)/(l - q)'^. 

From the above expression, it is easy to see that the coefficient a is arbitrary and the value of the index n depends 
on B. Taking qo = —0.81 ± 0.14 and jo ~ 2.16^o;7g[20] (the subscript denotes the present value), we can give the 
present range of B is 0.03 < Bq < 0.5. By calculations and analysis, the results of the expression (37) are as follows: 
when the real solution exists, the range of B is either B < ^'^^'^'^ or B > iI+|i/Z5 Considering 0.03 < Bo < 0.5, we 
find the range of B is 0.03 < B < ^'^'g^ and the index n are 7.36577 < n+ < 30.716 and 3.61737 < n_ < 5.26214. 
However, when there is not any real solution, the range of B is 11=1^ < B < 0.5 and the index n can be taken as 
any real number. 

We can point out that for the model of fi{R) = R + aR", /2 = 1 and G{L,n) = L,n, the corresponding results to 
the WEC given by us are just the same as the ones in Rcf. [16]. 

4. THE INSTABILITY OF GENERALIZED f{R) GRAVITY MODELS WITH ARBITRARY 

MATTER-GEOMETRY COUPLING 

Modified gravity must be stable at the classical and quantum level. There are in principle several kinds of 
instabilities to consider[21]. Dolgov-Kawasaki instability [22] is one of them. Below, wc will focus on this instability 
and generalize to f{R) gravity models with arbitrary matter-geometry coupling. 
The trace of the field equation (2) is 



(36b) 



(36c) 



(36d) 
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^ + 7F5g7j{2[/i - if[ + 2Gf;,)R] - 3nf[ - eDG/^ - AhiG'L„, - G)} 

57577 /2G'r, 



(38) 

where T = g^'^T^i,. As usual, we take fi{R) as fi{R) = R + e(p{R), where e must be small to compatibility with Solar 
System experiment [23]. Following[22], we expand the space-time quantities of interest as the sum of a background 
with constant curvature and a small perturbation: R = Rq + Ri, T = Tq + Ti, L = Lq + Li, and the space-time 
metric can locally be approximated by g^^ = rj^i, + hf^^^, where 77^11/ is the Minkowski metric. In fact, this is a local 
expansion over small space-time regions that are locally flat. Accordingly, /i (R) = R0 + R1+ e(p(i?o) + ff'{Ro)Ri + 
f[ [R] = 1 + e(^'(i?o) + £v"{Ro)Ri + ... and the linearized version of the trace equation (38) in the perturbations yields 

[6G(L„)/^'(i?o) + 3eip"iRo)]Ri - 3eip"{Ro)V'Ri + [12LoG'(L™)/^'(i?o) 
+12LiG'(L™)/^'(i?o) + 12G(L,„)i?o/r(^o) + 6eRo^"'{Ro)]Ri 
-6e(/?"'(i?o)Vi?i • Vi?o - [1 - 2G(L„0/^(i?o) + eip'{Ro) - Roe^"iRo) 

-moe^"'{Ro) + 3e^"'(i?o)V2i?o]i?i = 6/^(i?o)V2G(i™) (39) 
-/2(i?o)(4Li - Ti)G'(L„0 - 6(Lo + i'i)/^(i?o)G'(i,„) - [6G(L™)/^'(i?o) 
+3e(^"(i?o)]i?o - 12(4 + Li)i?oG'(L„0/^'(i?o) - 12Loii/^(i?o)G"(i™) 
+3e(^"(i?o)V2i?o, 

where V, and ovcrdot denote the gradient, Laplacian operators in Euclidean three-dimensional space and differ- 
entiation with respect to time, respectively, and the zero order equation 

/2(i?o)ToG'(L™) = -4/2(i?o)G(L,„) - 2e</,(i?o) + 4G(L„0i?o/2(^o) ^^^^ 
+4/2(i?o)ioG'(L„0 + 2i?oe(^'(i?o) - i?o - 2G(L™)i?o/^(i?o) - Roeif'iRo) 
has been used. By further calculation, the effective mass TOc// of the dynamical degree of freedom i?i can be given as 

mlf^ = [6G(L,„)/^'(i?o) + 3e^"(i?o)]-M2G(L„0./^(i?o) - 1 - eip'{Ro) ^^^^ 
+Roe^"iRo) + 3i?oe</'"'(i?o) - 3e^"' iRo)V^ Ro]- 

The dominant term on the right hand side is [6G(L,„)/2'(i?o) + 3e(/3"(i?o)]^^ and the effective mass squared must be 
non-negative for stability. Therefore, /{'(i?) + 2G(L,„)/2'(i?) > is the stability criterion for the generalized f{R) 
gravity models with arbitrary matter-geometry coupling against Dolgov-Kawasaki instabilities. 

Note that by taking G{Lm) = L,n and rescaling the function f2{R) as 1 + Xf2{R), this criterion is the Dolgov- 
Kawasaki criterion in f{R) gravity with non-minimal coupling between matter and geometry, which is just the same 
as the one in Ref. [24]. While by setting f2{R) = 1 and G{Lm) = ^m, the results given by us is the Dolgov-Kawasaki 
criterion in f{R) gravity with non-coupling, which is consistent with the results given in Ref. [25]. 

For models given by Eq.(27), the Dolgov-Kawasaki criterion is 

en{n - l)i?" + 2G(L,„)am(TO - l)i?™ > (42) 

where 



(-l)"e, ifi?<0 
e, if i? > 



(-l)"a, if i? < 
a, if i? > 



(43) 
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It is clear that the stabihty eriterion of these models don't relate to the values of j and s, and the space-time only 
depends on R. When n — m the inequality (42) gives, e + 2G(L,„)a > 0. 

It is worth stressing that the inequality (42) also can be obtained from the inequality (28) by taking 

A^^ = 0, (44a) 
B"^ = 0, (44b) 

^ = nin -l)a\R 1™ G'iL„,). (44c) 

From the above discussions, we find that for models given by Eq.(27), the energy conditions, the Dolgov-Kawasaki 
criterion and the condition for attractive gravity have the the same type of inequalities, but note that they are 
independent each other. 

From the inequality (42), we see that the viability of the model with respect to the Dolgov-Kawasaki instability 
criterion will depend not only on the value of the constants e and a, but also on the space-time metric under 
consideration. This fact will give further constraints on the Ricci scalar. 

5. SUMMARY 

So far, we have derived the energy conditions (SEC, NEC, DEC, WEC) in the generalized f{R) gravity models 
with arbitrary coupling between matter and geometry. For the SEC and the NEC, the Raychaudhuri equation, which 
is the physical origin of them, has been used. From the derivation, we found that equivalent results can be obtained by 
taking the transformations p — > p'^^^ and p p'^^^ into p + 3p > and p + p> 0. Thus by extending this approach to 
p — p>0 and p > 0, the DEC and the WEC in the generalized f{R) gravity models with arbitrary coupling between 
matter and geometry can be obtained. The condition to keep gravity attractive and the Dolgov-Kawasaki criterion 
in the generalized f{R) gravity models have been also given, but the approaches of deriving them are different. 

It is worth noting that the energy conditions and the Dolgov-Kawasaki criterion obtained in this paper are quite 
general, which include the corresponding results given in Refs.[15, 16, 24, 25] as well as in the general relativity (GR) 
as special cases. 

Furthermore, in order to get some insight on the meaning of these energy conditions and the Dolgov-Kawasaki 
criterion, we have applied them to a class of models. In these models the energy conditions, the Dolgov-Kawasaki 
criterion and the condition for attractive gravity have the the same type of inequalities. By analysis, we find that 
the Dolgov-Kawasaki instability criterion depends not only on the value of the constants e and a but also on the 
space-time metric under consideration. 

In addition, we have considered the special model with fi{R) = R, hiR) ~ cuR'^ and G{Lm) ~ L„i ~ —p. By 
virtue of the WEC and the present astrophysical observations, the values of parameters a and n can be constrained 
in this model. Of course, we will continue to study other models of f{R) gravity with arbitrary coupling between 
matter and geometry in our following investigations. 
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